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Abstract 

This paper presents a rational approach to modelling the triple velocity corre- 
lations that appear in the transport equations for the Reynolds stresses. All ex- 
isting models of these correlations have largely been formulated on phenomeno- 
logical grounds and are defective in one important aspect: they all neglect to 
allow for the dependence of these correlations on the local gradients of mean 
velocity. The mathematical necessity for this dependence will be demonstrated 
in the paper. The present contribution lies in the novel use of Group Represen- 
tation Theory to determine the most general tensorial form of these correlations 
in terms of all the second- and third-order tensor quantities that appear in the 
exact equations that govern their evolution. The requisite representation did 
not exist in the literature and therefore had to be developed specifically for this 
purpose by Professor G. F. Smith. The outcome of this work is a mathematical 
framework for the construction of algebraic, explicit, and rational models for 
the triple velocity correlations that are theoretically consistent and include all 
the correct dependencies. Previous models are reviewed, and all are shown to 
be an incomplete subset of this new representation, even to lowest order. 

I Introduction 

The state of the art in closure of the Reynolds-averaged form of the Navier-Stokes 
equations requires the solution of a differential transport equation for each nonzero 
component of the Reynolds stress tensor ( upu ] )• These equations are exact but contain 
a number of correlations which are unknown and must first be modelled in terms of 
known or knowable quantities. One such correlation, and the focus of this work, is 
the time-averaged triple velocity product ( uKFfUk ) whose spatial gradients represent 
the rate at which the Reynolds stresses are transported by the turbulent fluctuations. 
Conventional models of these correlations have been formulated with this physical 
role in mind. Specifically, and by analogy with the gradient transport hypothesis, 
Daly and Harlow (1970) assumed that the triple correlations are proportional to the 
spatial gradients of the Reynolds stresses to obtain 

s-, k dv~v~ 

-UiUjUk = C s -u k ui— 1 ) 

e dxi 
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where k is the turbulence kinetic energy (i.e., half the trace of ), e is its dissipation 
by viscous action, and C s is a proportionality coefficient. The Daly-Harlow model is 
the most widely used model in advanced simulations methods. However, while the 
triple velocity correlations themselves are independent of the order of the indices, the 
representation (the right-hand side of Equation 1) is not invariant under the permu- 
tation of the same indices. This is not usually a problem in thin shear layers where 
the flow direction is known a priori. It does, however, become so in general flows 
where turbulent mixing can be significant in all coordinate directions. Several alter- 
native models have been reported in the literature, and discussion of some of these is 
deferred until Section III. Suffice it to say here that these models, while being consid- 
erably more complex than Daly and Harlow, have not produced tangible or consistent 
improvements in the quality of predictions. The lack of improvement is equally true 
in both simple, two-dimensional free shear flows as well as in rapidly evolving flows 
in the presence of complex strain fields. This is the case, in our opinion, because all 
proposed models fail to reflect the dependence of the triple velocity correlations on 
an essential mechanism for turbulence transport, namely the spatial gradients of the 
mean velocity. That a rational model for Iquyufc should exhibit explicit dependence on 
the gradients of mean velocity will lie demonstrated conclusively in the next section. 
The provision of a theoretical framework for the formulation of such a model and a 
demonstration of how this explicit dependence can be introduced without recourse to 
empiricism provide the principal motivation for the present work. In Section II, the 
mathematical basis of the problem is presented. Our proposals for the formulation of 
a rational model for the triple velocity correlation are put forward in Section III and 
are discussed there in relation to existing models. Conclusions are given in Section IV. 
In developing a rational model for the triple velocity correlations, we sought to find 
a general representation of a third-order tensor in terms of both second- and third- 
order tensors. No such representation could be found in the published literature, and 
the need therefore arose for the development of the mathematical representation to 
specifically meet this requirement. This challenging task was performed by Professor 
G. F. Smith; his original work, which will no doubt find application in other fields of 
science, is included in the appendix. 


II The Mathematical Context 

We consider here the turbulent flow of a viscous, incompressible fluid with con- 
stant properties. The governing field equations are the Navier-Stokes and continuity 
equations, which are given by 

Dii; dp d 2 iii . . 

-a=-^ + "&r (2) 


duj 

dxi 


= 0 , 


( 3 ) 


where ii t is the instantaneous velocity vector, p is the static pressure and v is the 
kinematic viscosity of the fluid. In equations (2) and (3), the Einstein summation 
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convention applies to repeated indices, and D / Dt is the convective derivative defined 
as 

d _ d d 

Dt dt :> dxj 

The instantaneous velocity and pressure are decomposed into mean and fluctuat- 
ing parts as follows: 


(4) 


U t = Ui + Ui, p = P + p. 


(5) 


By substituting (5) into equation (2) and applying the Reynolds averaging rules (cf. 
Hinze 1975), there results the Reynolds-averaged Navier-Stokes equation: 


DUi 

Dt 


op d 2 Ui 

+ v- 


dx; 


dx] 


duiUj 

dxi 


( 6 ) 


Due to the nonlinearity of the convection term in equation (2), the time-averaging op- 
eration introduces additional unknowns in the form of the Reynolds stress tensor ujuj. 
This is the closure problem referred to in the Introduction. The most advanced meth- 
ods for resolving this problem rely on the solution of differential transport equations 
for UiUj. Derivation of these equations proceeds by multiplying the i - th component 
of equation (2) by uj and the j-tli component by u t . By summing and applying the 
Reynolds averaging rules, we have 


Convection 

D u i ll , 


Dt 


Production 


du 


dlh 


{ll II, jh + UjUk — ) 


dxi 


Dif fusion 


d_ 

dxi 


1 


du^j 


UiU,U k + -(p'lliSjk + p’lljSik ) - v 

P dx k 


Dissipation 

9 ( dui duj ' 

\ dx k dx k J 


Redistribution 


p 


+ - 


du i dui 

+ 


P \ dx i 


dx ; 


(7) 


The Reynolds averaging process has clearly resulted in additional unknown turbulence 
correlations, another manifestation of the closure problem. The literature contains 
numerous proposals for approximating these correlations (Speziale 1991 provides a 
comprehensive review of the subject). Our interest here is in the triple fluctuating- 
velocity components (run.,. The exact equation governing the evolution of UiUjUk 
reads (cf. Chou 1945): 


Du n a., 
Dt 


UiUjUl 


dlh 

dxi 


+ UjUkUl 


dh 

dxi 


+ U pupil 


djh 

dx i 
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II 


- UjUj- \-UjU k - \-U k Ui- 

\ dx k dxi ( 


_du k ui 


Jduiui 


+ UiUj— 1 - UjU k — 1 - U k Ui 

V dxi dxi 


dujm 
dx i 


d d 2 UiUjU k 

- — UiUjU k Ul + V — — 


( dui diij du k dui 

- 2 v U k -- 1 - Uj — — 

\ OX I OX I OX I OX I 


duj du k 
' dx i dxi 


It is immediately clear from line I of this equation that the gradients of mean velocity 
appear explicitly in the exact equation for the triple velocity correlations. Further, 
and as was shown by Chou f 1945), the gradients of mean velocity must also appear ex- 
plicitly in the correlation between the Reynolds stresses and the pressure fluctuations 
(line II above). Thus, by differentiating equation (2) to obtain a Poisson equation 
for the instantaneous pressure, and then by subtracting the mean, Chou obtained an 
equation for the fluctuating pressure: 

i dp d (u m u k Um^k) -y d l du k 

P d X k d X m dx k dx m 

Then, with the assumption of homogeneous turbulence, and for regions of the flow 
remote from solid boundaries, the following approximate solution was obtained: 


dp dp dp 

UiUj- b UjU k - b U k Ui—— 

dx k dxi oxj / 


+ Cijk 


The forms of the tensors, b nm ij k and Cij k , that appear in equation (10) are not material 
for the present purposes (they are given by equation (5.8) of the original reference). 
What is important here is the presence of the mean velocity gradients in the solu- 
tion. It therefore follows that a rational model for the triple velocity correlation must 
contain an explicit dependence on the gradients of mean velocity so that it remains 
consistent with the exact equation, and so that it stands a chance of representing the 
correct physics. Mathematically, this proposition implies that a model for the triple 
velocity correlations must be of the following functional form: 


i b ra n • H <i 


II i U j 11. k — F{jk 


where S mn and W mn are, respectively, the mean rate of strain and the mean vorticity 
tensors: 


Q — _ 

u mn ^ 
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= 1 fOUrn _ OUA 
2 Y 3x n 3x m J 


(13) 


The dissipation rate e is included in the argument of the general functional representa- 
tion in order to build up the turbulence length and time scales (/ and r, respectively) 
from the relations: 


l oc 


fc 3 / 2 


k 

T OC — 

e 


(14) 


Inspection of the Daly and Harlow (1970) model, equation (1), shows that it follows 
the incomplete representation 


,y ( duiUr 

UiUjUfo — -^ijk I 


(15) 


This is also true of all the models that have appeared in the literature. Hanjalic 
and Launder (1972), for example, derived a model from simplification of equation 
(8), but although they recognized that a complete model for II should contain the 
gradients of mean velocity, they proceeded to neglect this contribution from their 
workings. Thus, while their model represents an improvement over that of equation 
(1) (in being symmetric under the permutation of indices), its use in actual flow 
simulations has not produced tangible improvements. A discussion of some of the 
published proposals for the triple velocity correlations is included in the next section. 
Suffice it to say that all turn out to be of the functional form given by equation 
(15) rather than the mathematically consistent form of equation (11). In the next 
secion, we put forward our proposals for the formulation of a rational and consistent 
alternative. 


Ill The Present Proposal 

One way to proceed in the formulation of a rational model for EfuyuL is to attempt 
a term-by-term approximation of the various unknown turbulence correlations that 
appear in equation (8). The obvious difficulty with this approach lies in formulating 
adequate approximations to the large number of additional unknown terms that are 
involved. An example is the term that involves the gradients of the quadruple velocity 
correlation. An exact equation for this term is again not difficult to construct but will 
involve yet higher-order unknowns that will need to be approximated. Truncation of 
the process of introducing increasingly higher-order correlations can only be achieved 
by arbitrary considerations — something that is not in keeping with the present 
requirement for a rational formulation. 

Instead, the present proposal is to use Group Representation Theory to construct, 
in accordance with the functional relation of equation (11), the most general algebraic 
representation for the symmetric third-order tensor uiujuk in terms of both second- 
and third-order tensors. Both symmetric and antisymmetric second-order tensors 
are to be included in the representation. This representation is also required to be 
invariant under the three-dimensional full orthogonal group O 3 . 
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Surprisingly, the required mathematical representation does not exist in the lit- 
erature. Pennisi (1992) proposed a tensorial representation for a functional of a 
third-order tensor whose dependent variables consist of first- and second-order ten- 
sors. This representation, is clearly inadequate for the present purposes as it would 
lead to a model for the triple velocity correlations that does not contain dependence 
on the gradients of the Reynolds stresses. This dependence is required by term II in 
equation (11). 

The necessary representation therefore had to be developed specifically for this 
work. This task was achieved by Professor G F Smith. As it forms a distinct con- 
tribution to the subject, details of the theory are reproduced fully in the appendix. 
The complete representation for the triple velocity correlations that is consistent with 
the functional of equation (11) admits tensor dependencies up to fourth order. We 
shall henceforth dismiss terms that are of orders three and four because of the large 
numbers involved, equation (A4), and more importantly, because these terms do not 
introduce tensor parameters that are not already present in the linear and bilinear 
representation. 

With this simplification, the general representation of the triple velocity correla- 
tions suggested by Smith’s theory consists of the following group of terms: 

Group I: Terms that are linear in ( upi, ) tk . There are three linearly independent 
terms of this type, as indicated by (All). 

Group II: Terms that are bilinear in (u Ju]),k and muj. There are eleven such terms 
(A12). 

Group III: Terms that are bilinear in i u u and Sij. There are eleven linearly 
independent terms of this type (A13). 

Group IV: Terms that are bilinear in (u t Uj) : k and Wij. There are five such terms 
(A14). There are thus thirty terms in all, the combination of which produces a general 
representation of the form 

30 

u iUj u k = ]T (16) 

7=1 


where the (7 = 1,... ,30) are the third-order symmetric tensor-valued isotropic 
functions listed in A11-A14. The c 7 are material constants. 

We shall now consider the relation between this general representation and various 
ad-hoc models for the triple correlations. The literature contains a diverse assortment 
of such models, but these fall into two categories which are either linear or bilinear 
in the tensor dependencies shown in equation (15). 

First we consider models that are linear in the Reynolds stresses uluj. Examples 
of these models are the proposals of Mellor and Herring (1973) (henceforth referred 
to as MH) and Shir (1973). The MH model is given by 


k 2 

— UiUjUk = Cmh — 

e 


dujUk 

dx t 


+ 


du k u t 

dx 3 


+ 


duiUj \ 
dx k ) ’ 


(17) 


where Cmh is a, constant. Shir’s model contains only the third term in equation (17) 
and therefore does not exhibit correct symmetry. Both models are a sub-set of the 
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general representation in equation (16). The MH model, for example, corresponds 
to the X) Ajk term of equation (All). If all linear terms were to be included in a 
model, then two additional terms would appear, namely Ai pp 8j k and J2 A ppi 8 :) k . A 
consistent linear model for the triple velocity correlations would thus read 


UiZljUk — 


+ 

+ 


Mellor and Herring 




k 2 ( diijUk du k Ui duiUj 


dxi 


+ 


dx , 


+ 


dxi 


k 2 ( duiu p duju p 

Cmh. 2— I — — -8 lk + ,, Ski + 


c 


k 2 ( dk 


dx 


j'j’P % ! duku p 

On 


dx r 


MH, 3" 


^ dk_ 
dxP k ' dx , 


8ki + 


dk N 

dx ; , i r 


(18) 


Incidentally, this more complete formulation was also considered by Mellor and Her- 
ring (1973) who discarded it on the grounds of computational expediency. 

Next, we consider models that are bilinear in the tensor quantities. The Daly 
and Harlow (1970) model, equation (1), is of this category but, as already mentioned 
in the Introduction, is not independent of the order of indices and cannot therefore 
be included in the tensor representation shown in the appendix. 

A more complicated, but tensorially admissible expression has been proposed by 
Hanjalic and Launder (1972) (hereafter HL). They arrived at their model by com- 
bining lines III and IV of equation (8) after a proposition by Millionshtchikov (1941) 
and by neglecting lines I and II entirely. In addition, the transport term on the 
left-hand side of the same equation to obtain the following algebraic representation: 

k ( duju k du k Ui duiUj 

-UiUjU k = C hL~ UiUl —X f UjUl— b u k ui— 

t y dx; dx t dxi 



where Chl is a constant. This model corresponds to the Y^Aij p T p k term given in 
(A12). Note that the MH model can be viewed as an isotropized version of the HL 
model. 

Another model which is bilinear in the tensor quantities is that of Lumley (1978). 
He considered weakly anisotropic and inhomogeneous flows and applied stochastic 
considerations to obtain 

i j, r 

U t U j H;, — ( i | Qijk T f I, ( QirnrnAy k T Q mjm8ik T Gmmk8ij ) 

In terms of the notation used here, Gijk = J2 Aij p T pk , and the model given in equation 
(20) can be rewritten as 



, k 


( 21 ) 


Equation (21) shows that Lumley’s model is tensorially equivalent to that of Hanjalic 
and Launder in being bilinear in the tensor product of uiu m and duiu m /dx n . Both 
models, however, are only a subset of the 11 possible bilinear basis tensors (see (A12)). 
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Finally, Cormack et al. (1978) formulated a model in terms of the deviatoric. part 
of the Reynolds stress tensor, — 2/3 Sijk. In the mathematical framework outlined 

here, their model reads 


k / \ 

-UiUjUk = [C'CLS E] Aijk + C CLS E/ AppiSjk + CcLS-Aippfijk) 


+ 7C CLS E AippTjk + CcLS y ^ AppiTjk + CcLS E -AppqTq i 8 j k ^ 


(22) 


These six closure constants are linear combinations of the original four constants 
proposed by Cormack et al. (1978). Magnaudet (1993) also considered the Cormack 
model and provided a cross-diffusion extension which included the gradient of the 
isotropic dissipation rate. These terms were not included in the functional form 
proposed here, equation (11), as there is no evidence from equation (8) that such 
terms should be present. 

It should be recognized that since the tensor representations given in (All) - (A14) 
do not imply any asymptotic ordering of the terms, all these models should be viewed 
as arbitrary truncations of the full representation basis. Table I gives a summary of the 
models discussed and the corresponding basis terms involved. It is clear from Table 
I that the majority of these models are bilinear in the Reynolds stresses and their 
gradients. None of them contain the dependence on the mean velocity gradients that 
is required by the exact equation. It is not surprising, therefore, that there is little real 


Table I. Summary of u t UjUk models and corresponding tensor bases 



Model 

Basis 

Hanjalic and Launder 

Mellor and Herring 

Lumley 

Cormack et al. 

Magnaudet 

Tensor 

(1972) 

(1973) 

(1978) 

(1978) 

(1993) 

y ^ 


0.073 


0.069 


/ , A ipp ^ j k 




-0.136 


/.Appt&jk 




-0.632 


/ , ^ ipp T jk 




0.102 

0.16 

/ , Appi Tjf~ 




-0.068 


/ , AijpTpk 

0.11 


0.098 


0.125 

/ . AipqTpqS jh 



0.013 



/ j AppqTqiSjf- 




0.192 



difference in their performance across a wide range of benchmark turbulent shear flows 
(see, e.g., Schwarz and Bradshaw 1994, and Cormack et al. 1978). There are a number 
of flows in which the ratio of turbulence kinetic energy production to dissipation 
rates is near unity, and hence the gradients of the triple velocity correlations make 
only a modest contribution to the overall budget of the Reynolds stresses. In such 
flows, defects in modelling the triple velocity correlations are easily concealed by 
shortcomings in the closure of the other unknowns in equation (7). In other flows, such 
as in the self-similar wake behind a bluff body, for example, the ratio of production 
to dissipation is only about 10 percent, with the remainder of the balance provided 
by the gradients of the triple correlations. It is precisely in these “weakly sheared” 



flows that current turbulence closures seem to fail very badly. 

While the formulation of a triple-velocity correlations model that can be used in 
practical computations is outside the scope of this work, it is nevertheless useful to 
illustrate how such a model could be formulated. It is reasonable to argue that terms 
that are bilinear in (uiu]) t k and ujuj should be included in the compact model for the 
triple correlations, simply in order to ensure that these quantities will vanish as either 
the Reynolds stresses themselves, or as their spatial gradients, go to zero. Dependence 
upon the gradients of mean velocity can be introduced through the retention of a term 
that is bilinear in the Reynolds stress gradients and the mean rate of strain tensor. 
The resulting model would then read 


UiUjUk 



, c du k Ui t c 
“T ‘~P3 q x T Jpk 


du iu j 


dx,. 


(23) 


Alternatively, an argument can be made to the effect that the contributions to the 
triple-velocity correlations arising from .5'q and ITq are of the same order and, as 
such, can be combined. Thus, with the retention of a pair of terms that are bilinear 
in ( UiUj) } k and Sij or ITq there results 


UiUjUfz 


+ 



dujiik du k ui du{Uj 

Ujllp A llfcllp- 


dU p dujuk dU p du k iii 
dxj dXr 


+ 


dxp J 

dUp duiUj 
dxk dxr 


(24) 


The illustrative models represented bv equations (23) and (24) are, to our knowl- 
edge, the first rational, algebraic and explicit representations for the triple velocity 
correlations that are formulated from theoretical considerations and that allow for 
the explicit dependence of these correlations on the gradients of mean velocity. The 
coefficients cq and cq can be taken as constants and assigned numerical values by 
reference to the extensive data sets now available from Direct Numerical Simulations 
(e.g. Mansour et al. 1988). 


IY Closure 

The purpose of this paper was to advance a rational approach to modelling the 
triple velocity correlations of turbulence. A rational model was obtained by recourse 
to Group Representation Theory and by the development, specifically for the present 
purposes, of a representation for a symmetric third-order tensor in terms of both 
second- and third-order tensors, which are invariant under the full orthogonal group. 
Analysis of previous models showed all to be an incomplete subset of the full for- 
mulation that lacked, in particular, an explicit dependence on the gradients of mean 
velocity. This dependence appears explicitly in our representation and can there- 
fore provide the basis for the development of a rational model for the triple velocity 
correlation. A demonstration of how this model may be formulated was provided. 
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Appendix 


Oil the Representation of Third-Order Tensors 

G. F. Smith 

Department of Mechanical Engineering & Mechanics 
Lehigh University, Bethlehem, PA 18015, USA 


The following notation is adopted in this appendix: 


Tpq = UpUq 

Apq r = ( UpUq ) ^ 

Fijk — 7 / i li j U l z 

We consider the problem of determining the form of a third-order symmetric tensor- 
valued polynomial function 


Fijk •Fijk{Apq r . I pqi Spq,W pq) , ( A 1 ) 

which is invariant under the orthogonal group O 3 and which will be referred to as 
an isotropic function. The group O 3 is comprised of all 3x3 orthogonal matrices Q 
which satisfy 


QQ y = Q t Q = I = diag(l, 1 , 1), 

del Q = ±1. 

(32) 

The tensors appearing in (Al) satisfy 



Cijk — Cjki — Ckij — Cikj — Ckji 

Idjik , 

(.43) 

Aijk — Ajik, Tij — Tj i , S%j — S ji 5 

Wij = Wji. 



We may readily determine the number n 1 , 112, 113, 71,4 of linearly independent third- 
order tensor- valued functions which are of degree 1, 2, 3, 4 respectively in the compo- 
nents of the Aijk , . . . , Wij and which are invariant under O3. The procedure employed 
is similar to that adopted by Smith (1968). We have 

7i\ — 3, n 2 = 27, 723 = 225, 714 = 1429. (A4) 


The large number of terms of degree 3 and 4 in the Aijk , • • • , Wij leads us to limit 
consideration to the terms of degree < 2 appearing in (Al). Let 


B\) ] = 






rA 2 ) _ C.. 

D IJ ~ °IJ1 


B * 3) = W ■ 


(A5) 


The terms of degree < 2 appearing in the polynomial expression (Al) may then be 
written as 


c 


ijk 


— C'ijk T 


3 

\ ' 


/-)(«) rA) \ T? A 

U ijklm D lm > ^ijklmn^lmnA- 


Of — 1 


, 0) rO) n{P) 1 n[ a ) a r(«)_i_ 

' / j L ijklmnp^lm ^ np ' / > ^ ijklmnpq / *-i r mn J -'pq 1 

a, (3=1; a<.(3 


a = 1 


+ II, 


ijklmnpqr -^Imn-^pqr • 


(A6) 
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The tensors Cijk, D^l lm> . . . , H ljk i mnpqr appearing in (A6) are referred to as prop- 
erty tensors. The requirement that (A6) be invariant under O 3 imposes restrictions 
on the property tensors. Thus, the tensors Cijk, D^ lm , . . . must be invariant under 
0 3 , i.e., 

C ijk QipQjqQkrCpqr 

(-47) 

Dijklm = QipQjqQkrQlsQmtDpglst, (o! = (l, 2, 3),... 

must hold for all Q = ||Q; P |[ belonging to O3. Such tensors are referred to as isotropic 
tensors. There are no tensors of odd order which are invariant under O 3 . Even order 
tensors which are invariant under O 3 are expressible in terms of tensors of the form 
Sij, SijSki, Sij8kiS mn , ■ ■ • where Sij denotes the Kronecker delta tensor defined by 


8ij = 1 if ij = 11, 22, 33; 

S i:j = 0 if ij = 12, 21, 13, 31, 23, 32. 


(A8) 


See, for example, Smith (1994, p. 94) or Smith (1968). Thus, the tensors of odd 
order C ijk , D { °l lm , H i]klmnpqr appearing in (A6) are null tensors, i.e., all 

components of these tensors are zero. The expression (A6) then reduces to 


Cijk - EijklmnAlmn + % mnpq Al mn T pq + 

_l_ /^( 2 ) A c_L A 

' ^ijklrnnpq^lmn^pq ' ^ ijklmnpq^lmnW PQ * 


(.49) 


We may readily show that the number of linearly independent third-order symmetric 
tensor-valued isotropic functions which are linear in Aijk, bilinear in Aij k and T tJ , 
bilinear in A ljk and Sij, and bilinear in A ljk and W t] are given by 3, 11, 11 and 5 
respectively. We shall employ the notation E to indicate the sum of the three 

tensors obtained by cyclic permutation of the subscripts ijk on the summand. For 
example, 

E A t j k — Aij k T Ajki T A k ij , 

(410) 

E AijpTpk — AijpTpk T AjkpTpi T Aki p T p j . 

The three linearly independent terms appearing in (A9) which are linear in Aij k are 
given by 

El Aijk , E AippSjk-, y ^ AppiSji~. (411) 

The eleven linearly independent terms appearing in (A9) which are bilinear in A{ jk 
and Tij are given by 


T pp E A 


ijk , E AippTjk, AppiTjk, AijpTpk, E(A>j T Ajpi) T pk 


E A 


ipqT~pqSjk, ApgiTpqSjk, E A pp qTqiS j k , 


(412) 


^2 -AqppTqiSjki Tqq 22 AppjSjk^ Tqq 22 AjppSjk. 
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The eleven linearly independent terms appearing in (A9) which are bilinear in Aijk 
and Sij are given by 


$pp X^ A ijk 5 X^ Aipp Sj k 7 XZ App{S jk , X^ AijpSpk^ X^ T Ajpi ) 

X^ -A-ipqSpqS jk-) E A-pqiSpqS X^ AppqSqiSjk ? (A13) 

E A-qppS qi$ S ' qq X^ AppiSjk^ S qq X^ A{ppSj} 


The five linearly independent terms appearing in (A9) which are bilinear in and 
1 TX are given by 


E E \-A-ipj T Ajpi) E Aipq ^^pq$ jki 

E Appg Sj Jz ? E AqppSqiSjk. 


(.414) 


The polynomial expression (Al), which is invariant under Os and in which terms of 
degree three or greater in the A^k, . . . , W t j are neglected, is then given as a linear 
combination of the 30 terms listed above, i.e., 

30 

Cut = E chj! (.415) 

7—1 


where the (7 = 1, . . . , 30) are the third-order symmetric tensor-valued isotropic 
functions listed in (All),. . . , ( A 1 4 ) . The c 7 are material constants. 
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